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We show how to include in the existing calculations the corrections to the isovector coupling 
arising in chiral effective field theory recently found in Ref. 0. The dominant effect can be taken 
into account by conveniently redefining the static spin matrix elements (S p , n ): the largest one is 
reduced, on average, a 10%. To show the impact of these corrections we recalculate the limits on the 
WIMP-proton spin dependent cross scetion set by COUPP. We also give practical formulas to obtain 
(Sp, n ) given the structure functions in the various formalisms/notations existing in literature. We 
argue that the standard treatment of the spin-dependent cross section in terms of three independent 
isospin functions, 5*00(9), Sn(q), Soi(q), is redundant in the sense that the interference function 
5"oi (?) is the double product |Soi(g)| = 2y/S 00 {q) \/Sii(q). 



PACS numbers: 95.35.+d, 12.60.Jv 



I. INTRODUCTION 

The formalism for the cross section of the elastic scat- 
tering of weakly interacting massive particles (WIMP) 
off nuclei, was formulated many years ago for the neu- 
tralino [H Q , a well motivated supersymmetric dark mat- 
ter candidate Q. The dominant interactions of a non- 
relativistic WIMP with a nucleus are the coherent spin- 
independent and spin-spin interactions [U, [E[ • Nowdays 
many other particle physics models that predict a dark 
matter candidate have been proposed. However, in the 
non-relativistic limit, only a limited number of opera- 
tors remain, as has been shown with the help of model 
independent effective filed theory [(|. To calculate the 
contribution to the cross section of these operators, the 
corresponding nuclear responses must be evaluated. 

The scope of this paper is to clarify some previously 
unrecognized properties of the response functions of the 
spin dependent formalism and show how to take into ac- 
count the effects of low energy pion physics recently found 
in Ref. Q. 

In Section [XT] and Appendix [A] we show that the in- 
terference function Sqi is not an independent function 
but that it is determined by the isoscalar and isovector 
functions Sqo and Six- 

In Section Mil we show how to account for the strong 
interactions correction and also compare, with the aid 
of Appendix [B] the results of different calculations for 
nuclei of experimental interest. 

To illustrate the impact of these effects we recalculate 
the limits on the WIMP-proton cross section using the 
last data of the COUPP experiment ||. 



II. ONLY TWO FUNCTIONS 

The effective axial interaction of a spin 1/2 Majorana 
fcrmion with nucleons 



reduces in the non-relativistic limit to the operator 
(2V2G F )a N 4m x m N (4S x ■ S N ) and leads to a WIMP- 
nucleus cross section that, after Refs. [l], Hj], is given in 
the form 
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S(q) = a^Sooiq) 



S(q) 
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- a aiS 01 (q) + afSn(q). 



(1) 
(2) 



Although nuclear physics calculations are carried out 
in the isospin representation defined by the isoscalar cou- 
pling ao = a p + a n , and the isovector coupling a\ — 
o-p — a„, for practical purposes is natural to think in 
terms of neutrons and protons. The quantities usu- 
ally employed by experimentalists are thus the couplings 
a p = i a o + a i)/2, a n = ( a o — a i)/2, and the spin ma- 
trix elements of the protons and neutrons group in the 
nuclear ground state with total angular momentum J: 
(Sp) = (J, M = J\J2f=iSf\J,M = J), (Sn) = (J,M = 
J\ 53j=i Z I J, M = J). In the protons-neutrons repre- 
sentation, Eq. ([2]) at zero momentum transfer reduces to 



S(0)= 2J+1J+ T 1 \a p (S p )+a n (S n )\ 2 . (3) 

The differential cross section is usually then rewritten as 
da tua 



7 a A m SD (q), 



(4) 



dE R 2n\v 2 

where the total WIMP-nucleus cross sections at q = is 



<ta(0) = 8^^-5(0), 
and the form factor normalized to one at q = 



f> SD (q) 



s(oy 



(5) 



(6) 



The functional form of the structure functions Sij ob- 
tained from shell-model calculations is exactly a polyno- 
mial times an exponential in the dimensionless variable 
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y = (qb/2) 2 (or u — 2y), being b the oscillator size pa- 
rameter, if the single particle states used to build the 
nuclear wave function are wave functions of the three- 
dimensional harmonic oscillator potential. This is not 
the case if the eigenfunctions of the Wood-Saxon poten- 
tial are used [9h11|] and when the effects discussed in 
Section |Hl] are included; anyway, it is always possible 
to fit them with a polynomial or an exponential times a 
polynomial in y or u. Hence in general 
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or Si. 
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(7) 



It thus seems that three nuclear response function in 
the isospin representation Sn, Soo, ^01 are necessary to 
furnish the cross section at finite momentum transfer, 
while only the two numbers (S Pi7l ) are necessary at q = 0. 
This asymmetry is unnatural, q = is just a value of the 
variable q, and actually it does not hold. 

The functions Soa(q), S%i(q) are monotonic decreasing, 
always positive at every q, while Sqi is always positive 
or always negative at every q. The simple reason behind 
this behavior, as shown in Appendix [B] is that we can 
define two functions Sq (q) and Si (q) (that maintain their 
sign) such that the isoscalar and the isovector functions 
Sqo and Sn are nothing but the square of So and Si, 



Soo(g) = So(q), Su(q) = Sf(q), 



(8) 



and the interference function is the double product, 
Soi = 2S (q)Sx(q), \S i(q)\ = 2^ S 0Q (q)S 11 (q), (9) 

with sign poi = s 8' n (Soi) given by the relative sign be- 
tween So(q) and S±(q). 

The physics at q = is determined by the coefficients 
of order zero of the polynomials in Eq. (J7J, 

Soo(0) = 4°), Sn(0) = c£\ poi = sgn( C (° 1 ) ). (10) 

Comparing Eq. ([2]) evaluated at q — with Eq. ([3]) we 
have 



7T J 

2J+U+1 



(S oo (0) + Sii(0)±Soi(0)), 



being the plus sign associated with the protons and the 
minus sign with the neutrons. The static spin matrix 
elements, in the light of Eqs. dHJ) , ([9]) , (fT0|) , are then given 
by 



(V c oo +PoiV c ii )> 




where A is + or - in both, depending on the given nucleus. 



III. INCLUDING STRONG INTERACTION 
CORRECTIONS 

In Ref. it has been shown that strong interaction 
effects due to pion exchange that arise in chiral pertur- 
bation theory renormalize the isovector coupling a\. The 
authors also present a new shell-model calculation for the 
Xenon isotopes 129 Xe, 131 Xe giving numerical fits of the 
functions Sij that include these new corrections. 

Decomposing the spin-spin operator as a sum of the 
longitudinal and transverse electric operators projections 
of the axial current, see Appendix [B] it is shown in 
Ref. 0| that the isovector coupling ai is renormalized as 
5j-a\ in T Gl 5 where 



<y r = i-2^- + <yoi, 



and Scai in C 5 with 

^^TxpnF^q 



Sr = 1 - 



2c 3 pq 2 



2m N g A (ml + q 2 ) F 2 (4ml + q 2 ) 



(12) 



+ *o 1 . (13) 



The first q 2 -dependent term of Sj- and 8c arise at the one- 
body level (lb) currents, while the second q 2 dependent 
term in 5c and the momentum independent term 



5a\ = — ^rI(—k-\ ' 

F£ 3 GrriN ' 



(14) 



is due to two-body (2b) currents. Here — 92.4 MeV 
is the pion decay constant, = 138.04 MeV the pion 
mass, rriN the nucleon mass, g npn = 13.05 the pion- 
nucleon coupling, A = 1040 MeV a cut-off scale, p the 
nuclear density, p £ [0.10,0.12] fm~ 3 , I ~ 0.58 - 0.60 is 
a dimcnsionless factor. The largest uncertainties are in 
the parameters C3 € [—2.2, —4.78] and k — (2c3 — C4) £ 
[7.2,14.0], see Ref. [lj]. Considering the possible range 
in p, I and fc, and calling for simplicity 5a± = 5 we have 
for the interval and the average value: 

5 e [-0.135,-0.314], and <5 av = -0.224. (15) 

In principle, the calculation should be redone for all the 
nuclei of experimental interest. However we now show 
how to include them in existing calculations. The argu- 
ment is the following. 

The effect of the 2-body currents at q = can be trans- 
ferred into the static spin matrix elements by replacing 
a\ in S(0) with ai(l + 5). Using the defining relations of 
do,i we find 

S(0) 2 b = 2J + 1 J + T 1 |Qp(S p ) 2 b + an( S n)2b| 2 , 



with 



(S p ) 2b = {S p ) + -«S p > - <S B », 

(S„>2h = <S B )-|«S p )-(S„)). (16) 



3 



Obviously the total contribution of the proton and neu- 
tron groups to the total angular momentum does not 
change, (S p ) 2 b + (S„)2b = (S p ) + (S„). Note that also 
the sign of the spin matrix elements is important in the 
determination of the size of the correction. 

In Table U we report the spin values for the most im- 
portant isotopes employed as detecting medium in ac- 
tual experiments and their values after the inclusion of 
the 2-body current effects using the average value S av . 
The starting values are taken from papers that claim to 
present the best shell-model calculation for the given iso- 
tope (we rounded to 3 decimal digits where necessary). 

The case of 19 F needs a particular comment. The ex- 
periments employing 19 F, COUPP 0, SIMPLE [11 and 
PICASSO [13] , employ the spin values of Ref. jig, of- 
ten quoting the compilation of Ref. [16], (S p ) = 0.441, 
(S n ) = —0.109. The successive refined shell-model cal- 
culation of Ref. .17] using the more realistic Wildhental 
interaction found (S p ) = 0.4751 and (S„) = -0.0087. 
The protons contribution is thus similar but the neutrons 
contribution is now a factor ~ 50 smaller and not 4 times 
smaller than the protons one. 

The new calculations of Ref. Q allows us to clarify the 
question. These authors use the same Wildhental inter- 
action, the harmonic oscillator wave functions and the 
full s — d shell structure, thus the results should agree 
with the ones of Ref. [17J. They furnish the analytical 
expressions of the response functions for 19 F but do not 
calculate the static spin values. We use Eq. (|B3|) to ob- 
tain Sii(0) and Eqs. ([TO]) . (fTTjl to get the spin values. We 
find (S p ) = 0.475(5) and (S„) = -0.008(68), in perfect 
numerical agreement with the values of Ref. [l7J. Such 
an agreement is also found in the case of 23 Na, see Fig. 0] 
in Appendix IB] 

In Ref. [6j new shell-model calculations for 73 Ge, 127 I, 
129,131 Xe are also given. We have verified that in these 
cases the spin values differ largely from the other known 
calculations, see also Ref. [l9j]. Given that the same au- 
thors in Ref. [f| state explicitly that these calculations 
are "exploratory", we do not use in Table HI 

We see from Table U that with 5 av the dominant spin 
matrix element (that determine if the isotope is sensitive 
to the spin-dependent WIMP-proton or WIMP-neutron 
cross section) typically suffer a 10% reduction while the 
value of the smaller spin matrix element increases but 
always remain more than one order of magnitude smaller 
than the dominant one. 

At finite momentum transfer there are two energy 
scales to be compared. Experiments that measure the 
recoil energy typically restrict the search window below 
100 keV to reduce the background. Using the oscillator 
size parameter [f| 0, H3 



b = 



41.467 



45v4- 1 /3 - 25A- 2 / 3 



fm, 



100 keV corresponds to y ~ 0.07 in the case of 19 F, y ~ 
0.39 for 73 Ge and y ~ 0.83 in the case of 131 Xe. The other 
scale is set by the momentum transfer of the order of the 



Table I. Spin matrix elements for the odd-mass isotopes em- 
ployed in dark matter search experiments. The values of the 
last two columns include the 2-body currents corrections as 
given in Eq. with 8 = <5 a v of Eq. (fT5|). 


Isotope (J) [Ref.] 


(S P ) 


(Sn) 


(S p )2b 


(S,j)2b 


19 F (1/2) [17] 


0.475 


-0.009 


0.421 


0.045 


23 Na (3/2) [17] 


0.248 


0.020 


0.222 


0.045 


29 Si (1/2) [17] 


-0.002 


0.133 


0.013 


0.118 


73 Ge (9/2) [10] 


0.030 


0.378 


0.069 


0.339 


127 I (5/2) [11] 


0.309 


0.075 


0.283 


0.101 


129 Xe (1/2) [7] 


0.010 


0.329 


0.046 


0.293 


131 Xe (3/2) [7J 


-0.009 


-0.272 


-0.038 


-0.242 


133 Cs (7/2) [18] 


-0.318 


0.021 


-0.280 


-0.017 



pion mass, 100 MeV, that corresponds to y ~ 0.2, 0.28 
and 0.38, respectively, q 2 dependent corrections can thus 
be neglected for light nuclei but should be considered in 
medium-heavy and heavy nuclei. 
We anyway note that 



S£(0) S£(0) 



(17) 



as shown in Fig. [T] in the case of 131 Xe, blue full and 
dashed lines. The correction 2 /A 2 is practically negligible 
while the lb correction in aj- is already included in the 
standard calculations 0, H in the form 1 / (m 2 + q 2 ) using 
the Goldberger-Treiman relation g^NN — 9A"mN / F v . 

The normalized isospin functions have the property, as 
discussed in details in Refs. I20L I21L 



$oo 0?) = 



■Soo(g) 
5*00 (0) 



*u(g) 



Sn(O)" 



(18) 



For light nuclei they are practical identical for all the val- 
ues of interest for the experiments, in heavy nuclei the dif- 
ference get larger increasing q 2 . This fact is consequence 
of the isospin symmetry that treats proton and neutrons 
on an equal footing. At small y, once normalized the 
isospin form factors behave like 1 — ay + by 2 + 0(y 3 ). In 
light nuclei as 19 F, 23 Na, 29 Si, the polynomials are of low 
order, typically up to fourth degree, thus this relations is 
respected on a larger range of y, that anyway cover the 
range of momentum transfer caused by dark matter elas- 
tic scattering. For medium heavy and heavy nuclei the 
more complicated structure requires higher order polyno- 
mials, from sixth to nine-th degree, thus the differences 
are more accentuated. 

To conclude, the corrections can be taken into account 
in the standard expression of the cross section, Eq. (Q}, 
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Figure 1. Isospin form factors for the isotope 131 Xe built with 
the functions of Ref. See Eq. (|18jl for definitions. 



using the corrected static spin values (S Pi „) 2 b in 0^(0), 
Eq. ([5]), and using the known Su function in the form 
factor © normalizing them to one by factoring out their 
zero momentum transfer value. 

In Fig. [5] we show the impact of the 2b currents on 
the limits of COUPP (CF 3 I). We use the standard 
Maxwellian distribution truncated at v esc — 544 km/s, 
the velocity of the Sun vq = 230 km/s, the velocity of 
the Earth v E = 244 km/s and p = 0.3 GeV/cm 3 for 
the local dark matter density. The efficiencies, exposures 
and a threshold are specified in Ref. [8|. The two red 
curves correspond to two different models for the bub- 
ble nucleation efficiency of fluorine. Note that we find 
the agreement with the published curves using an upper 
limit on the number of events of N UL = 14.5, that is the 
90% confidence level limit obtained with the Feldman- 
Cousins method with 13 observed events (7 of 20 nuclear 
recoil events removed by the time isolation cut) and an 
expected background of 4.5. 

The black dashed lines are the new limits calculated 
with <5 av . The reduction in the mass region above the 
minimum is about 10%. The effects should be included 
also in analysis concerning the impact of nuclear physics 
uncertainties on the reconstruction of the WIMP proper- 
ties, mass and cross section, from the experiments with 
spin-dependent sensitivity [22j. 




10 

nfy. (GeV) 

Figure 2. The solid red lines are the COUPP @] 90% con- 
fidence level corresponding to two different models for the 
bubble nucleation efficiency of fluorine 0. See the text for 
further details. The black dashed line are calculated with 
corrected spin values of Table [T] 



Appendix A: The isoscalar-isovector interference 
function 



The spin-spin operator 



A 1 

O(q)=4S x .^-(a l 

»=i 



(Al) 



with Si and the spin and coordinates of the z-th nu- 
cleon Ta\p) — \p), Tz\n) = —\n) and 1, the identity op- 
erator in isospin space, is expressed in the standard for- 
malism [l[ (also in Ref. 0) as the sum of the transverse 
electric and longitudinal projections of the axial current 

O(q)=4S x -(T el5 (q)e ± +£ 5 (q)e ). 

eb is a spherical unit vector in the direction of the quan- 
tization axis taken to be one of q and e± are the unit 
vector orthogonal to this direction. These operators are 
decomposed into multipoles as 



i=l 



■ ai 5rff)T(q), 



A 1 

^ 5 = E^( 



o l + aiS c ff)L(q), 



(A2) 



(A3) 
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where the detailed expressions for T(q) and L(q) are un- 
necessary for our argument. The coefficients Sj- and 5c 
account for the effects discussed in Section Hm Given the 
nuclear wave function \A), then 

\M\ 2 ex (\(A\\T &l5 (q)U)\ 2 + \(A\\£ 5 (q)\\A)f). (A4) 

There is no T el 5 — £ 5 interference in the modulus squared 
because of the orthogonality but in Eq. (|A4[) each term 
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Figure 3. The isoscalar-isovector interference function Sen 
and the double product of Eq. ([9]l. Functions for 73 Ge are 
from Ref. Hp, for 127 I from Ref. [ll[ (Bonn A set), for 131 Xe 
from Ref. [7| (lb set). 



Figure 4. The structure functions for 23 Na. The superscript 
V and H indicate the that the functions are built from the 
functions of Ref. [l7| and of Ref. 0] respectively, using the 
prescriptions of Appendix [B] The superscript RD indicate 
the polynomial expressions of Ref. [llj . 



in Eq. (|A4I) gives contributions proportional to ag, a\ 
and the interference a^a\. The functions Sij(q) arise as 
a sum of pieces after a rearrangement of these terms: 

\M\ 2 oc (a 2 S Q0 (q) + afSniq) + a axS i(q)). (AS) 

On the other hand, equivalently, the spin-spin operator 
can be decomposed, as done in Refs. |17| . 

0(q) = 4S X ■ [a S + aiSi), 



where 

1 A 

S o(?) = -^lS,e"^ r % (A6) 

i=l 
1 A 

i=l 

are the nuclear spin operators associated with the isospin 
couplings. Performing the multipole decomposition on 
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these operators we are led to 

\M\ 2 cx \a (A\\S (q)\\A) +a 1 S 1 (A\\S 1 (q)\\A)\ 2 . (A7) 

Now all the corrections are included in S\ and we can 
define two functions such that Eq. (|A7|) takes the form 

\M\ 2 cx (a S 2 (q) + ai S 2 (q) + 2a a 1 S (q)S 1 (q)). (A8) 

Eq. (|A5|) and Eq. (|A8|) are equivalent: the isoscalar and 
the isovector functions Sqo and 5*11 are squared of some 
functions So(q) and Si(q) and the interference function 
is the double product 2So(q)Si(q). 

We explicitly show this fact in Fig. [3] for three iso- 
topes used in actual experiments. The functions refer- 
ring to 73 Ge are taken from Ref. [13], the ones for 127 I 
from Ref. [ll[ (Bonn A potential set) and for 131 Xe from 
Ref. 0- See also Fig. HI bottom panel, for 23 Na. The 
small differences that one can find in some cases can be 
ascribed to the intrinsic uncertainties of the numerical 
fits and numerical rounding. 



Appendix B: Three equivalent notations 

In the notation of Ref. [13] we can identify So = £Iq, 
Si = 17 1. From these the normalized form factors are 
defined, F i3 -(g) = f2?,(g)/fi?,(0). 



In the notation of Ref. |6| the basic spin response func- 
tions are, N,N' — p,n, F%,' N (corresponding to the ax- 
ial transverse electric operator T cl 5 ) and F^,', N (corre- 
sponding to the axial longitudinal operator C 5 ), thus 

F£ N ' = ±(F»' N '+F»,r). (Bl) 

The subscript "44" designate the spin-spin interaction. 
From these, the form factors in the isospin representation 
are 

^44 - 7^44 +^44 + ^44 )■> 

F\\ = \(F!f + Fit - 2Fit). (B2) 

The standard isoscalar and isovector functions Sqo and 
Sn read: 

2 7+ 1 2 7-1-1 

S u (q) = —J—n 2 (0)F u (q) ee A^-FH(q). (B3) 

l07T 7T 

This equivalence is illustrated in Fig. |3] for the iso- 
topes 23 Na, that to our knowledge is the only nucleus 
for which the calculation in the three formalisms exists 
in literature. The results of Ref. [13] and Ref. [3] coin- 
cide exactly for both Soo and Sn. The isoscalar function 
of Ref. [n| agrees with the others two, while it is not 
the case for the isovector function. Note that these last 
ones are furnished as polynomial fits that are valid up to 
y ~ 0.35, for y > 0.4 they become negative. The double 
product rule is respected by the three calculations. 
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